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Abstract 

We propose a modified Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 
for the summation of large ln(l/x), x being the Bjorken variable, which con- 
tains an extra dependence on momentum transfer Q compared to the con- 
ventional BFKL equation. This Q dependence comes from the phase space 
constraint for soft real gluon emissions in the derivation of the BFKL kernel. 
The modified equation gives a slower rise of the gluon distribution function 
with 1/x for smaller Q, such that the predictions for the structure func- 
tions F2^l{x, Q"^) involved in deep inelastic scattering match the HERA data. 
These results are similar to those obtained from the Ciafaloni-Catani-Fiorani- 
Marchesini equation, which embodies both the InQ and ln(l/x) summations. 
The reason that the conventional BFKL equation can not be applied to the 
low Q (soft pomeron dominant) region is supplied. 
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I. INTRODUCTION 



The deep inelastic scattering (DIS) experiments performed at HERA [|T| 
have stimulated many intensive discussions on the behavior of the structure 
function -^2(0:, Q"^) at small Bjorken variable x for various momentum transfer 
Q. One of the important subjects is to explain the data of F2 using the 
known evolution equations. In the low x region with ln(l/x) ^ InQ, the 
appropriate tool is the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 
which sums large ln(l/x) in the gluon distribution function to all orders. 
Since the BFKL equation does not depend on momentum transfer explicitly, 
its predictions for the gluon distribution function and for F2 exhibit a weak 
Q dependence. However, the HERA data show that the rise of F2 at small x 
varies with Q in a more sensitive way: The ascent is rapid for large Q, which 
is attributed to hard (or BFKL) pomeron contributions. The slow ascent at 
small Q, which can not be explained by the BFKL equation, is interpreted 
as the consequence of soft pomeron exchanges. 

To understand the Q-dependent rise of the gluon distribution function 
and the structure function, the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equation for the InQ summation should be included some- 
how. One may resort to the DGLAP equation directly, because the relevant 
splitting function Pgg oc 1/x also gives an increase at small x. However, 
the DGLAP rise is milder, and a steep nonperturbative input of the gluon 
distribution function at low Q must be adopted. An alternative is the pt- 
factorization theorem |Q, which has been shown to be equivalent to the 
collinear (mass) factorization, the framework on which the DGLAP equa- 
tion is based. By including the next-to- leading ln(l/a;) from singlet quark 
contributions, the rise becomes faster, and a flat input is proposed 0. The 
double scaling tendency exhibited by the predictions from the DGLAP equa- 
tion with the next-to-leading logarithms taken into account and by the data 
of F2 has been explored 0. The conclusion also favors a flat nonperturbative 
input. The above investigations imply that the structure function at low Q 
is dominated by soft pomeron contributions, and the large-Q behavior is the 
result of the DGLAP evolution. 

Another possibility is the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) 
equation [0 , which embodies the BFKL equation for small x (and intermedi- 
ate Q) and the DGLAP equation for large Q (and intermediate x). To derive 
the CCFM equation, the angular ordering of rung gluons in a ladder diagram 
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is assumed, which can be regared as a mixture of the rapidity ordering for 
the BFKL equation and the transverse momentum ordering for the DGLAP 
equation. The CCFM equation has been studied in details in . It is found 
that the resultant gluon distribution function shows the desired Q depen- 
dence, and the predictions for F2 are in agreement with the HERA data. 
This equation, though considering only the gluon contributions, predicts a 
steeper rise than the DGLAP equation does, since the running coupling con- 
stant involved in the splitting function is evaluated at the gluon transverse 
momentum, instead of at 

The above analyses, though consistent with the data, did not address the 
issue why the BFKL equation can not be applied to DIS in the low Q region, 
where it is supposed to be more appropriate because of ln(l/x) ^ InQ. In 
this paper we shall derive the BFKL equation in a more careful way from 
the viewpoint of the Collins-Soper-Sterman resummation technique [^. By 
truncating the loop momenta of real gluons at Q in the evaluation of the 
evolution kernel, a modified equation is obtained. Since real gluon emissions 
are responsible for the BFKL rise, the rise is rendered slower at smaller 
Q, where the phase space for real gluon emissions is more restricted. We 
argue that this is the reason the conventional BFKL equation, overestimating 
real gluon contributions, fails in the low Q region. The modified equation 
is simpler than the CCFM equation, and explains the HERA data of the 
structure functions F2^l well. 

In Sect. II we derive the Q-dependent BFKL equation using the re- 
summation technique, and show that it approaches the conventional one as 
Q ^ 00. In Sect. HI the evolution of the gluon distribution function in 
X predicted by the modified BFKL equation is compared with those by the 
DGLAP equation, the conventional BFKL equation, and the CCFM equa- 
tion. The modified BFKL equation is solved numerically, and the results of 
the gluon distribution function and of the structure functions F2^l are pre- 
sented in Sect. IV. Section V is the conclusioin, and the Appendix contains 
the details of the involved calculations. 



II. FORMALISM 



In pH] , 0] we have proposed a unified derivation of the known evolu- 
tion equations that sums various large logarithms based on the resummation 
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technique [p|, such as the DGLAP, BFKL, CCFM, and Gribov-Levin-Ryskin 



(GLR) equation |12|. The GLR equation includes the annihilation effect of 
two gluons into one gluon in the region with Q and 1/x being simultaneously 
large. Briefly speaking, the resummation technique relates the derivative of a 
parton distribution function to a new function involving a new vertex |13 



The new function is then expressed as a factorization formula involving the 
subdiagram containing the new vertex and the original distribution function. 
Evaluating the subdiagram according to the specific orderings of radiative 
gluons simply, the derivative of the distribution function becomes the corre- 
sponding evolution equation. The complexity of the conventional derivation 
of the evolution equations is thus greatly reduced. 

In this section we apply the above formalism to DIS of a proton with 
light-like momentum = P^5^^ in the small x limit, where x = —(^j {2P ■ 
q) = Q'^/{2P ■ q) is the Bjorken variable, q being the momentum of the 
photon. The unintegrated gluon distribution function F{x,kT), describing 
the probability of a gluon carrying a longitudinal momentum fraction x and 
transverse momentum pt, is defined by 

x{P\F^iy-,yT)F^^O)\P) , (1) 



in the axial gauge n ■ A = = 0, n = being a vector on the light cone. 
The ket \P) denotes the incoming proton, and is the field tensor. Averages 
over the spin and color are understood. To implement the resummation 
technique, we allow n to vary away from the light cone (n^ 7^ 0) first. It 
will be shown that the BFKL kernel turns out to be n-independent. After 
deriving the evolution equation, n is brought back to the light cone, and the 
definition of the gluon distribution function coincides with the standard one. 
That is, the arbitrary vector n appears only at the intermediate stage of the 
derivation and as an auxiliary tool of the resummation tecchnique. 

Though F at X ^ does not depend on P"*" explicitly as shown later, it 
varies with P"*" through the momentum fraction implicitly, which is propor- 
tional to {P^)~^. Hence, the derivative P^dF/dP^ can be written as 

.F{x,Pt)^-x^F{x,Pt). (2) 



dP+ ' ' ' dx 
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Because of the scale invariance of F in n, as indicated by the gluon propa- 
gator, {-i/P)N'"'{l), with 



= r - ^ + n'r^, , (3) 



F must depend on P via the ratio {P-nf' / n^. We have the chain rule relating 
P+d/dP+ to d/dn 

Vf3 = Sj3+ being a vector along P. The operator d/dnp appplying to the gluon 
propagator gives 



d 



—irN^'^' + r'N''^). (5) 



dnp n ■ I 

The loop momentum l'^ [l'^') contracts with a vertex in F, which is then 
replaced by a special vertex 

f/3 = 7 • 6) 

V ■ nn ■ I 

This special vertex can be read off the combination of Eqs. @) and (^. 



The contraction of I" hints the application of the Ward identities |Tn, 111 



^ ^1 1.'). ~ K?. /K , n-?. ' y' ) 



for the quark-gluon vertex, and 



F {k + lf 



(8) 



for the triple-gluon vertex T ^^x. Similar identities for other types of vertices 
can be derived easily. Summing all the diagrams with different differentiated 
gluons, those embedding the special vertices cancel by pairs, leaving the one 
where the special vertex moves to the outer end of the parton line We 
obtain the formula, 

-xA-F{x,pt) = 2F{x,pt) , (9) 
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described by Fig. 1(a), where the new function F contains one special vertex 
represented by a square. The coefficient 2 comes from the equahty of the 
new functions with the special vertex on either of the two parton lines. Note 
that Eq. is an exact consequence of the Ward identity without approx- 
imation 0. An approximation will be introduced, when F is related to F 
by factorizing out the subdiagram containing the special vertex, such that 
Eq. (H) reduces to a differential equation of F. 

It is known that factorization holds in the leading regions. The lead- 
ing regions of the loop momentum / flowing through the special vertex are 
soft and hard, since the factor 1/n ■ I with 7^ in Eq. (^) suppresses 
collinear divergences 0. For soft and hard /, the subdiagram containing 
the special vertex is factorized according to Figs. 1(b) and 1(c), respectively. 
Fig. 1(b) collects the soft divergences of the subdiagram by eikonalizing the 
gluon propagator, which will be explained below. We extract the color factor 
from the relation fabcfbdc = —NcSad, the indices a,b, . . . being indicated in 
Fig. 1(b), and Nc = 3 the number of colors. The corresponding factorization 
formula is then written as 



F,{x,pt) = tN^g'J —--T.^xVpHN''^ 



2m5{l')F{x + Ipt + ItI) + ^^f^F(x,pT) 

(10) 

where iN^. comes from the product of the overall coefficient —i in Eq. (P) and 
the color factor —N^ extracted above, g is the coupling constant, and the 
triple-gluon vertex for vanishing / is given by 



-g^^yxPx - guxxP^, + 2gxf,xP^ . (11) 



The denominator —2xP-l comes from the eikonal approximation {xP — l)"^ ~ 
—2xP-l. The ffist term in the brackets corresponds to the real gluon emission, 
where F{x + l^/P^, \pt + 1t|) implies that the parton coming out of the 
proton carries the momentum components xP^ + and + 1^ in order 
to radiate a real gluon of momentum /. The second term corresponds to the 
virtual gluon emission, where the 9 function sets the upper bound of It to 
Pt to ensure a soft momentum flow. 



6 



It can be easily shown that the contraction of P with a vertex in the quark 
box diagram the partons attach, or with a vertex in the gluon distribution 
function, gives a contribution down by a power 1/ s, s = {P + qY, compared 
to the contribution from the contraction with -0^. Following this observation, 
Eq. (|10D is reexpressed as 



dH 



vl 



T ''T_ 

P 



F{x,pt) 



27rz5(/2)F(x + /+/P+,|pT + lT|) 

(12) 



The eikonal vertex v^, comes from the term xP,y (divided by xP^) in Eq. (11 
and the eikonal propagator 1/v ■ I from l/{xP ■ I), which is represented by 
a double line in Fig. 1(b). The remaining metric tensor g'^'^ has been ab- 
sorbed into F. We apply the strong rapidity ordering x + /P^ ^ x to the 
evaluation of Eq. ([T2|), that is, approximate F{x + l^/P~^, \pt + 1t|) by its 
dominant value F{x, \pt + lr|)- The integration over and to infinity 



gives 



Fs{x,kT) = y 



—n 



\l+=oo 
\l+=0 



X 



TT 2n-[n+/^ + 2n-/+2] 
F(x, \kT + h\) - 0{kT - It)F{x, kr) 

T 



[^(^' |kT + ItI) - ^(4 - lT)Fix, kr)] , (13) 



with as = Ncas/TT and n = {n~^,n~,0) has been assumed for convenience. 
The first line of the above formulas demonstrates explicitly how the n de- 
pendence cancels in the evaluation of Fig. 1(b). 

It can be shown that the contribution from the first diagram of Fig. 1(c) 



vanishes like 11 





r d% 


' 1 


2 J 


TT 


Jjrp 


1 


xP'^u 



In 



ixP+vY - xP+v 



(14) 
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at X — >■ 0, with the factor z/ = y (f ■ n)^/n^. This is the reason F does not 
acquire an exphcit dependence on the large scale P"*" as stated before, and 
the transverse degrees of freedom must be taken into account, leading to the 
PT-factorization theorem ||^. The n dependence residing in Fig. 1(c) then 
disappears with the vanishing of Eq. ([T^. If neglecting the contribution 
from Fig. 1(c), i.e., adopting F = F^, Eq. (|^) becomes 



d\n{l/x) 



^siPT) I ^ \f{x, Ipt + ItI) - e{pl - 11)F{x,pt)\ , (15) 



which is exactly the BFKL equation. The argument of has been set to the 
natural scale pt- It is then understood that the subdiagram containing the 
special vertex plays the role of the evolution kernel. Obviously, the BFKL 
kernel is gauge invariant after deriving the evolution equation. Now we make 
the vector n approach S'^', and the definition of F returns to the standard 
one. 

A careful examination reveals that it is not proper to extend the loop 
momentum of a real gluon to infinity when deriving Eq. (|13]), since the 
behavior of F{x + l^ / P^), vanishing at the large momentum fraction, should 
introduce an upper bound of To obtain a more reasonable evolution 
kernel, we propose to truncate /"*" at some scale, and a plausible choice of this 
scale is of order Q. The cutoff of introduces a n dependence. However, 
this n dependence can be absorbed into the cutoff, which will be regarded 
as a fitting parameter and determined by the data of F2 (for some value of 
Q). Then the gauge dependence does not appear explicitly in the modified 
evolution equation. It can be verified that the predictions are insensitive to 
the cutoff, and a good choice of the cutoff is Qji^v), where v is the gauge 
factor appearing in Eq. (11^). Performing the integration over l" and 
Eq. (O) becomes 



F,{x,pt,Q) 



a, r d'^lr —1 



d+=Q/{2u) 



x'F{x,\pT + \T\,Q)-e{p^T-^l)Hx,PT,Q)] , (16) 



where the Q dependence of F from the constraint of phase space for the real 
gluon emission has been indicated. With this modification, the real gluon 
has smaller phase space at lower Q. 
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We derive the modified BFKL equation 



^siPT) I '^[F{x,\vT + \TlQ)-e{pl-ll)F{x,PT,Q)] 
J TlLj. 

-<Pt)J- , (17) 

where the first and last terms correspond to the lower and upper bounds of 
respectively. Obviously, Eq. ( |T7|) approaches the conventional equation 
( p!5D in the Q ^ oo limit. The first term on the right-hand side of Eq. (p!7| ) is 
responsible for the rise of F at small x, and the last term acts to moderate the 
rise We then expect that the ascent of F is slower at smaller Q, for which 
the effect of the last term is stronger, and deviates from that predicted by 
the conventional BFKL equation. Note that the Q dependence is attributed 
to the constraint of the phase space for radiative corrections, instead of to 
the In Q summation. 



dF{x,pT,Q) 
ciln(l/x) 



III. EVOLUTION IN x 

Before proceeding with the numerical analysis, we extract the behavior of 
the gluon distribution function at small x from the modified BFKL equation, 
and compare it with those from the DGLAP, BFKL, and CCFM equations. It 
will be observed that the modified BFKL equation and the CCFM equation 
give similar evolution in x. It is most convenient to consider the Fourier 
transformed version of the equations. We reexpress Eq. (|1^) as |11 



dF{x,pT,Q) _ f cPl 
= a 



d\n{l/x) 



J -:^[F{x, Ipt + ItI, Q) - 0{Ql - ll)F{x,PT, Q)] 

'""'J- fTo^ • ^^^^ 



with a fixed coupling constant Og- Note that the loop momentum in the 
virtual gluon emission is truncated at the scale Qq, instead of at pt- This 
modification is acceptable, since the virtual gluon contribution plays the role 
of a soft regulator for the real gluon emission, and setting the cutoff to Qo 
serves the same purpose. Furthermore, the replacement of px by Qq allows 
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us to solve Eq. (|I8]) analytically, and the solution of F maintains the essential 
BFKL features. 

The Fourier transform of Eq. (]TB|) to the conjugate h space leads to 



with 

S{h, Q) = 2a,[ln(go&) + 7 - In 2 + K^i^Qh)] . (20) 

where the Bessel function Kq comes from the last term on the right-hand 
side of Eq. (p!8[), and 7 is the Euler constant. Equation ( |T9| ) is trivially solved 
to give 

6, Q) oc exp[-S(6, Q) ln(xo/x)] , (21) 

Xq being the initial momentum fraction below which the gluon distribution 
function begins to evolve according to the BFKL equation Transforming 
Eq. (pTD back to momentum space, we have 

F{x,pT,Q) = / bdbJo{pTb)F{x,b,Q) . (22) 
Jo 

Assuming that the above integral is dominated by the contribution from the 
small b region, F is expected to grow as 

F(x, 0, Q) oc exp[2a, ln(xo/x) ln(g/go)] . (23) 

This power- law rise with 1/x is steeper compared to the DGLAP evolution, 



exp 



A^cj^lj^ MQ/Aqcd) 



\ /3o X ln((5o/AQCD) 



(24) 



being the first coefficient of the QCD beta function, and is more sensitive 
to the variation of Q compared to the conventional BFKL evolution, 

exp[41n2a,ln(a;o/a;)] . (25) 

Remind that the DGLAP evolution with x is not strong enough, while the 
conventional BFKL rise is Q-independent. Equations (p3l)-(p5|) then make 
clear why the modified BFKL equation can explain the data well. 
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At the same time, we shall study the CCFM equation, whose Q depen- 
dence comes from the In Q summation. It is written as 

F{x,pT,Q) = F^^\x,pT,Q)+ C dz j ^e{Q-zq)As{Q,zq) 

Jx J TTO^ 



xP{z,q,pT)F{x/z,\YiT + {I - z)(ilq) , (26) 



with the splitting function 

P = as{pT 



h ^Ns{.z,q,PT)- + z{l - z] 

1 — z z 



(27) 



The so-called "Sudakov" exponential As and the "non-Sudakov" exponen- 
tials A are given by 



\s{Q,zq) 
^Ns{z,q,PT) 



exp 
exp 



-a. 



dp^ f^-PT/p dz' 



^0 dz' i-Pt dp^ 

Z' J{z'qy p"^ 



1 - Z' 



(28) 



The upper bound zq of the variable z' in Eq. (^) takes the values 

1 if 1 < (pr/q) 

zo = Pr/q if z < {pr/q) < 1 
z if {pt/ q) <z . 



(29) 



The Sudakov exponential A5 collects the contributions from the rung gluons 
obeying the strong angular ordering, and is the result of the In Q summation. 
Those gluons which do not obey the angular ordering are grouped into Aats. 
For the simpler derivation of the CCFM equation using the resummation 
technique, refer to In the small x region we adopt the approximate 

form of the CCFM equation [R|, 



F(x,pT,Q) = F'^'\x,pT,Q) + as{pT) f - f ^e{Q - zq)d{q ~ f,) 

Jx z J Tcq'^ 

xANs{z,q,pT)F{x/z, Ipt + (1 - z)q\,q) . (30) 



The extra theta function 6{q — fi) introduces an infrared cutoff of the variable 
q. The flat nonperturbative driving term F*^"^ will be defined later. 
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We extract the approximate CCFM evolution in x as a comparision. Ap- 
plying d / d\n.{l / x) to both sides of Eq. (|30D gives 

dF{x,pT,Q) _ dz r d^q 



as I — I ^e{Q - zq)e{q- ^)ANs{z,q,pT) 
Jx z J Tiq 



dhi{l/x) Jx z J nq 

Again, we assume a fixed coupling constant cxg- The derivative of the fiat 
driving term F*^") has been dropped. The differentiation of the lower bound 
of z vanishes, because of = 1) = 0. The arguments |pt + (1 — z)q\ and 
g of F in the integrand have been replaced by |pt + q| and Q, respectively, 
for simplicity. The identity dF/d\n{l/x) = dF/dlnz and the neglect of the 
derivatives of other factors in the integrand lead to 

dF{x,pT,Q) _ J d r d^q 



d\xi(l/x) Jx dz J nq"^ 

xANs{z,q,PT)F{x/z, \pt + q|,Q) 



j2 

a, I -\e{Q - q)e{q - |pt + q|, Q) • (32) 



To derive the last line of the above formulas, we have employed F{^ = 1) = 
and Ajvs(l, q',pr) = 1 according to Eq. (pH]). The Fourier transform of the 
above equation is given by 



dF 



dhi{l/x) 



2a, HQ/^i)F (33) 



in the small h region. It is obvious that the CCFM equation predicts a similar 
evolution in x to Eq. (^). This similarity will be verified numerically in the 
next section. However, its Q dependence is due to the transverse momentum 
ordering as indicated by the function 6{Q — q), a source different from the 
constraint on the longitudinal momenta of real gluons in the modified BFKL 
equation. 

Another interesting observation is that the CCFM equation gives a steeper 
rise at small x than the DGLAP equation does. The reason is attributed to 
the different choices of the argument of the running coupling constant in 
the splitting function, which is Pt in the former (see Eq. (pTj)) and q in the 
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latter. Hence, asipr) does not run in fact, as the variable q is integrated 
over in the CCFM equation (|30|) . While as{q) evolves to as{Q) when solving 
the DGLAP equation. Because of agipr) > c^siQ), the CCFM evolution is 
stronger. 



IV. NUMERICAL RESULTS 

We solve both the modified and conventional BFKL equations numerically 
following the repetation method proposed in by starting with a "flat" 
gluon distribution function p|, 



F^Hx,Pt) = ^N.il-xfexpi-p'^/Ql) , (34) 

with Qo = 1 GeV. The active flavor number n/ in the running coupling 
constant as is set to 4. The normalization constant A^^ will be determined by 
the data of the structure function F2(x, Q^) at some Q, and then employed 
to make predictions for other Q. We solve Eqs. (|15D and (0) with F^'^^ 
in Eq. (|3^ ) inserted into their right-hand sides. Substitute the obtained 
solution of F into the right-hand sides of Eqs. (|15D and (|17D, and solve them 



again. Repeat this procedure until a limit solution of F is reached. In the 
numerical analysis we introduce a lower bound p-r > 1 GeV to avoid the 
infrared divergence in agipr) from the pt diffusion at small x. The CCFM 
equation is also solved in a similar way with the driving term F^'^^ 

3 /.I 
F^'\x,pT,Q) = —,N,eM-pllQl) -eiQ-zpT) 

wich reduces to F^''-\ when the functions 6 and I^ns are set to unity. 

The expressions of the structure functions F2^Li according to the pt- 
factorization theorem 0], are given by 

_ 



F2(x,Q2) = / A r _PLH,[xli,pT.Q)F{i,PT,Q) + -FL{x,Q'), 

Jx ^ Jo TC 2 

^^^^ 

Fl{x,Q') = f'f r''^HL{x/^,pT,Q)F{^,pT,Q), (37) 

Jx E Jo TT 
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Pc and p[. being the upper bounds of px which will be specified later. The 
hard scattering subamplitudes if2,L denote the contributions from the quark 
box diagrams in Fig. 2, where both the incoming photon and gluon are off 
shell by — — and p^ — —p^, respectively, p = (^P"*", 0, Pt) being the 
parton momentum and ^ the momentum fraction. They are written as 



Pt 



z^ + {l- zY - 2z{l -2z)^ + 2z 



I , l + A 
A 1-A 



(38) 



and Hl = Hl+^ + HI, with 



-2e'^Czl--- + z^ — 
^""^n^'U z^'q^A^ 



2^ 
' A^ 



(39) 



the sum of the contributions from Figs. 2(a) and 2(b), and 

z 

2zx — z — 2z 



HI = 



-2et^Cz\ 

^ 71 



1 - 



p¥ 



2pI b'' 



{z — zx) 



g2 



A 



A I- A 



2x\l-2z 



1 



2.§) - .A^ 



A^ 



A^ l-A^ 



pI\1- 2^(3 - z)pI/Q^ + Qz^p^tIQ 



A^ 



X 



A" 



1 , l + A 

— in 7 

A I- A 



(40) 
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from Fig. 2(c). The factors A, B and C are A = J I - Az'^pH Q"^ , B 



1 — z — zp^/Q"^ and C = 1/2. The detailed derivation of iJ2,L is presented 
in the Appendix. We have assumed a vanishing charm quark mass for sim- 
phcity. To require if2,L be meaningful, the upper bounds of pt are set to 

Pc = min ^Q, , p'^ = min ^Q, ^-^Q^ ■ (41) 

It is easy to observe that the terms in the braces of H2 approaches the 
splitting function 

P,,iz) = C[z' + (1 - zf] (42) 

in the Pt ^ limit. 

We evaluate Eq. (^6]) for = 15 GeV^ with F{C,,pt,Q) obtained from 
the modified BFKL equation, and then determine the normalization constant 
as Ng = 3.29 from the data fitting. When Q"^ varies, we adjust Ng such that 
xg has a fixed normalization Jq xgdx. F2 for Q"^ = 8.5, 12, 25, and 50 
GeV^ are then computed, and results along with the HERA data [|l| are 
displayed in Fig. 3. It is obvious that our predictions agree with the data 
well. The curves have a slower rise at smaller Q, which is the consequence 
of the phase space constraint for real gluon emissions. For comparision, we 
show the results from the conventional BFKL equation, and from the CCFM 
equation. The former are almost Q-independent, and close to those from the 
modified BFKL equation only at large = 50 GeV^. Therefore, their match 
with the data is not very satisfactory, especially in the low Q region. The 
latter, also consistent with the data, are similar to those from the modified 
BFKL equation as expected. 

The normalization constants Ng for F2 are then employed to evaluate the 
structure function Fl{x,Q^) using Eq. (^) for the corresponding Q. The 
predictions for = 8.5, 12, 15, 25, and 50 GeV^ are presented in Fig. 4, 
which behave in a similar way to F2{x,Q'^), but are smaller in magnitude. 
Note that the uncertainty of the experimental data [|I^ shown in Fig. 4 is 
still large: The vertical lines adjacent to x = 10~^ for = 8.5 and 12 GeV^ 
represent the error bars. The central values slightly above 0.5 are outside of 
the plot range. Hence, it is not difficult to explain the current data. To draw 
a concrete conclusion, more precise data are necessary. 

The gluon density xg is defined by the integral of the unintegrated gluon 
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distribution function over pt [§, 



xg{x,Q^)= ( -^F{x,pT,Q). (43) 

Jo TT 

The dependnece of xg on x for = 8.5, 15, and 50 GeV^ derived from the 
modified BFKL equation is displayed in Fig. 5. If letting the parameter A 
characterize the behavior of xg ~ x~^ at small x, we find the values A ~ 0.3 
for = 8.5 GeV^ and A ^ 0.4 for Q'^ = 50 GeV^. The former is consistent 



with that obtained from a phenomenological fit to the HERA data |16]. The 
latter is close to that from solving the conventional BFKL equation numer- 
ically ll^, i.e., from hard pomeron contributions. As explained before, the 



smaller A at lower Q is due to the more restricted phase space for real gluon 
emissions. The fiat gluon density corresponding to soft pomeron exchanges 
can then be understood following this vein. 



V. CONCLUSION 

In this paper we have studied in details the modified BFKL equation 
proposed from the viewpoint of the resummation technique. This equation 
possesses a Q dependence from the constraint of the phase space for soft 
real gluon emissions. It gives rise to the desired behaviors in x and in Q of 
the gluon distribution function and the structure functions -F2,l, which are 
consistent with the HERA data. The predictions are close to those from the 
CCFM equation, whose Q dependence is, however, due to the InQ summa- 
tion. The rise of the gluon density at small x predicted by the modified BFKL 
equation is steeper than that by the DGLAP equation, and is more sensitive 
to the variation of Q compared to that by the conventional BFKL equation. 
We have also explained why the conventional BFKL equation overestimates 
the real gluon contributions and thus the rise of the structure functions in 
the low Q region. Real gluons, responsible for the BFKL rise, in fact should 
have more restricted phase space at smaller Q. 

This work demonstrates the simplicity and success of the resummation 
technique in the study of the DIS structure functions at small x. However, 
the power-law rise of the structure functions predicted by either the con- 
ventional or modified BFKL equation violates the unitarity bound |T^. A 
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further application of our formalism to the issue of the unitarity of the BFKL 
evolution will be published elsewhere |jl9[ . 



This work was supported by the National Science Council of Republic of 
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Appendix 

In this Appendix we extract the hard scattering subamplitudes from the 
box diagrams in Fig. 2. The photon and gluon with the off-shell momenta 

= —Q^ and p"^ = —px, respectively, annihilate into a quark with mo- 
mentum k and an antiquark with k'. These diagrams can be regarded as a 
parton-level cross section, which are then convoluted with the unintegrated 
gluon distribution function F as shown in Eq. (0). The structure functions 
-^2,L are extracted from the DIS hadronic tensor W^i, through the contrac- 
tions 

where the longitudinal structure function is, in terms of the standard ones 
Fi^2, written as F^ = F2 — 2xFi. 

In the infinite-momentum frame of the proton, the only nonvanishing 
component of P is P~^, and the gluon carries the parton momentum C,P~^+Pt- 
In the center-of-mass frame of the photon and gluon, the gluon momentum 
is transformed into p with the space component p in the z direction. The 
relation q + p = leads to 

q',=q'+h\' = -Q' + \p\'=pl-p'-Q'. (45) 



iFrom the parton-level Bjorken variable 

X Q2 



^ '^p-q 2(pogo + |p||q|) 



(46) 
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with Eq. ( ^51) inserted, we solve for po, which is expressed as 

Q 1 + 2zpyQ^ 



Po 



2 Jzil -z + zpyq^ 



(47) 



It is then easy to obtain 



Q ^i + 4zv/g2 
2" ^z{i - z + zpyq^) 

Q l-2z 



2 ^z{i -z + zpyq^) 

Since each of the quarks shares half of the total energy, we have 
, 1/ ^ , Q l-z + zp^Q' 

2 2 _ ^ + zpyq^) 



(48) 
(49) 



(50) 



Similarly, the invariants s = {p + q)"^, t = (jp — k'Y and u = {p — k)"^ are 
written as 



u 



■^(l + ,/l+422£_cos^^ 
2^ 



1 - \ l + Az'^^co^e 



(51) 



where 9 is the angle between the momenta k and p. 

We first compute —g'^^W^u, where W^i, is the parton-level hadronic ten- 
sor. Contracting —g^'^ to the box diagrams, the squared amplitude from 
Fig. 2(a) is given by 



{k - qf 



7m /^7^ 



(k-q) 



r7 



(52) 



where the factor 1/4 denotes the average over the spins of the photon and 
gluon, C = 1/2 is the color factor from the identity Tr{T°-T^) = (1/2)5'*'', 
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j^a,b iQQiYLg the color matrices, and is the charge of the quark q in unit of 
e. A straightforward calculation leads to 



(53) 



Since Fig. 2(b) is idential to Fig. 2(a) under the interchange of k and k', we 
have {—g^''Wfj,u)b = (— 5''^'^^/iz/)a(^ ^ u). A similar evaluation of Fig. 2(c) 
gives 



The complete expression is then written as 



tu 



(54) 



SiregagC 



t 



u 



tu 



(55) 



Next we calculate the contraction P'^P^Wfj^y. Note that the proton mo- 
mentum P is not parallel to the gluon momentum p, because the gluon is off 
shell, while the proton is on shell. In the center-of-mass frame P possesses 
different components from those in the infinite-momentum frame. From the 
invariants, 



P-P^PoPo-PsIpI = 0, 



2P-q 2{Poqo + Ps\q.\) ' 
the components Pq and P3 in the center-of-mass frame are solved to be 

Q 1 



(56) 



'^^^z{l-z + zpyQ^) ' 

Q 1 + 2 v/g^ 

2^ ^z{i -z + zpyQ^)^i + Az-'pyq^ 



(57) 



Because of P^ = 0, it is easy to obtain the transverse components 



Pt 



1 



(58) 
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It will be found that all the contributions from Fig. 2 to P'^P'^W^y can 
be expressed in terms of the invariant, 



P ■ k = Pq^o ~ -Ps^o cos 6 — Prko sin 6 cos ( 



(59) 



with 



Poko 
Prko 



Ax ' 

i + 2V/g^ 

4a; ,/l + 4zV/Q^ 



Jl + 42V/Q^ 



(60) 



being the azimuthal angle of the momentum k. Following a similar proce- 
dure, we derive the squared amplitudes from Figs. 2(a), 2(b) and 2(c), 



,2 Q' 



P^P^'W, 



P^P'^W, 



2{P ■ ky - —P ■ k 



87relasC-\2{Q'+p'){P-ky 



tu 

z 



2z 



2z 



[l-Acos9)+p'^ 



P-k 



(61) 



respectively, with A = ^Jl + Az'^p'^/Q'^. The complete expression is given by 
P'^P'^W^u = {P^P'W^,)^ + (P^'P'^W^,)^ + (P^P^#,,,)^ . (62) 



To extract the hard scattering subamplitudes, we integrate Eqs. ( |55| ) and 
( |6lD over phase space. For —g^'^Wp^p which is independnet of the azimuthal 
angle 0, the integral is written as 

H2 = / d^k'd^k5\p + q-k-k') 
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= ^ J dcosOzi-g'-^W,,) . (63) 

Substituting Eq. (|55|) into Eq. (^) and integrating it over 0086*, we obtain 
Eq. (^), where the relation = — has been inserted. For another con- 
traction P^P^W^u, which depends on 0, the integral is given by 

HL = -^j d cos ed(l)2z—P^P''W^, . (64) 
After a tedious calculation, we obtain Eqs. (^) and (^Of ). 
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Figure Captions 



Fig. 1. (a) The derivative —xdF/dx in axial gauge, (b) The soft structure 
and (c) the ultraviolet structure of the 0{as) subdiagram containing the 
special vertex. 

Fig. 2. The quark box diagrams contributing to the hard scattering sub- 
amplitudes H2^L- 

Fig. 3. The dependence of F2 on x for = 3,5^ 12, 15, 25, and 50 GeV^ 
derived from the modified BFKL equation, the conventional BFKL equation, 
and the CCFM equation. The HERA data are also shown. 

Fig. 4. The dependence of Fl on x for Q'^ = 8.5, 12, 15, 25, and 50 GeV^ 
derived from the modified BFKL equation, the conventional BFKL equation, 
and the CCFM equation. The HERA data |5| are also shown. 



Fig. 5. The dependence of xg on x derived from the modified BFKL 
equation. 
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